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QCD perturbation theory ignores the compact nature of SU( 3) gauge group that gives rise to 
the periodic 0-vacuum of the theory. We propose to modify the gluon propagator to reconcile 
perturbation theory with the anomalous Ward identities for the topological current in the 9- vacuum. 

As a result, the gluon couples to the Veneziano ghost describing the tunneling transitions between 
different Chern-Simons sectors of the vacuum; we call the emerging gluon dressed by ghost loops a 
“glost”. We evaluate the glost propagator and find that it has the form G(p) = (p 2 + Xto P /p 2 ) _1 
where \to P is the Yang-Mills topological susceptibility related to the v[ mass by Witten-Veneziano 
relation; this propagator describes confinement of gluons at distances ~ Xtop 4 — 1 &n. The same 
functional form of the propagator was originally proposed by Gribov as a solution to the gauge 
copies problem that plagues perturbation theory. The resulting running coupling coincides with the 
perturbative one at p 2 xtop , but in the infrared region either freezes (in pure Yang-Mills theory) 

or vanishes (in full QCD with light quarks), in accord with experimental evidence. Our scenario 
makes explicit the connection between confinement and topology of the QCD vacuum; we discuss 
the implications for spin physics, high energy scattering, and the physics of quark-gluon plasma. 

PACS numbers: 25.75.Bh, 13.87.Fh, 12.38.Mh 


QCD possesses a compact gauge group that allows for topologically non-trivial gauge field configurations. These 
configurations realize homotopy maps from the gauge group to the space-time manifold. For example, the homotopy 
map from the SU(2) subgroup of the gauge group to the Euclidean space-time sphere S 3 describes the instanton 
solution PP. However the compactness of the gauge group is ignored in perturbation theory, and this may be at the 
origin of problems marring the perturbative approach. 

In QCD, one of these problems is the existence of Gribov copies [2] - multiple solutions of the gauge-fixing condition 
that make the perturbative approach ambiguous. In Coulomb gauge, the emergence of Gribov copies can be traced 
back to the existence of energy-degenerate vacua with different Chern-Simons numbers [3]. A natural question arises 
- is it possible to formulate QCD perturbation theory in a way that is consistent with the topological structure of 
the theory? In this Letter we argue that the answer to this question is positive. We find that the resulting gluon 
propagator naturally describes confinement, i.e. non-propagation of color degrees of freedom, and the running coupling 
displays the screening of color charge at large distances. 

In Minkowski space-time, instanton solutions represent the tunneling events connecting the degenerate vacuum 
states with different Chern-Simons numbers 


X(t) = J d 3 x K 0 {x,t), (1) 

where K 0 is the temporal component of topological current 

K, = (d?A°' a + ± g C abc A p b A^ ; (2) 

the first term in I\q is the density of Abelian “magnetic helicity” while the second term is its non-Abelian generaliza¬ 
tion. 

The chiral anomaly in QCD leads to non-conservation of the axial current 

dp. Ja = 2N f Q(x) + ^(2im / )9 / 75Q/, 

/ 


( 3 ) 
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where to/ are the masses of quarks, Nf is the number of flavors, and 

Q(X) = ^ F ^ X)P ^ V( - X) (4) 

is the density of topological charge normalized by f d 4 x Q(x ) = is; for finite action field configurations is is an integer. 
The density of topological charge can be represented as a divergence Q(x) = cC K lt of the gauge-dependent current 

§• 

Veneziano [1] has demonstrated that the periodic 0-vacuum stucture in QCD can be captured by introducing a 
massless “ghost” in the correlation function of the gauge-dependent topological current ([2]): 

K^(q) = i J d 4 x e iqx (K»{x)K v { 0)) (5) 

where /x 4 = Xtop is the topological susceptibility of pure Yang-Mills theory. Note that the r.h.s. of Eq. ([ 5 ]) has the 
“wrong” sign, i.e. the ghost does not describe a propagating degree of freedom. This means that the ghost cannot be 
produced in a physical process; however the couplings of the ghost (that describe the effect of topological fluctuations) 
certainly can affect physical amplitudes. A similar “dipole” ghost had been earlier introduced by Kogut and Susskind 
[5] in the analysis of axial anomaly in the Schwinger model. This procedure has been demonstrated to solve the 17^(1) 
problem in QCD m- 


The physical meaning of Eq. ([5| becomes apparent if 
coordinate x(t) in a crystal [8]: 


one compares it to the correlation function of the electron’s 


i [ dt e iut (T{x{t)x{ 0)}) - 

J uA m 


1 d 2 E{k) 
u} 2 dk 2 


k =0 


( 6 ) 


where E = k 2 /2m* is the energy of an electron with an effective mass to* and quasi-momentum k in a crystal. The 
emergence of the pole in Eq. Q signals the possibility of electron’s propagation in the periodic potential of the crystal 
due to tunneling. Note that the pole emerges not just from a single tunneling event (corresponding to the instanton 
in QCD), but sums up the effect of many tunnelings throughout the crystalline lattice. 


The analogy between Eq. ([ 5 ]) and Eq. ([6]) can be made even more apparent if we choose the frame with q ^ = (w, 0) 
and use the (gauge-invariant) analog of coordinate given by Eq. ([l]). The expression Eq. ([ 5 ]) then takes the form 
completely analogous to Eq. <©: 


i 


dt e iu,t (T{A(t)A(0)}} 


«-*> _J_ d 2 E{0) 

oj 2 u) 2 d6 2 


(7) 


where V is the volume of the system, and E{6) = e(0)V is the energy of the Yang-Mills vacuum. The energy density 
of the vacuum e(0) is a periodic function of the 0 angle that is analogous to the quasi-momentum k in Eq. ([6]). At 
small 0, we can expand e(0) and write 


e( 6 ')=Ai 4 y, (8) 

which exhibits the physical meaning of /x 4 as of the topological susceptibility Xtop = /x 4 of the Yang-Mills theory; note 
that a term linear in 0 is forbidden by P and CP invariances of QCD. 


It is well known that topological susceptibility vanishes, order by order, in perturbation theory. Perturbative 
description thus corresponds to /x — > 0 in Eq. ([ 5 ]), or to the limit of the infinitely heavy electron, to* — > 00 in Eq. ©. 
Infinitely heavy electrons do not respond to electromagnetic fields, as the corresponding coupling of electromagnetic 
current j = e q/m* to the gauge field jA ~ 1/to* vanishes in the limit to* —» 00 . In this case the dynamics of photons 
is not sensitive to the periodic structure of the crystal, and one can build the usual perturbation theory of photons. 
However when m* is finite, and is of the order of the frequency of the external gauge field, photons can be absorbed 
and re-emitted by the electrons, and these processes severely affect the photon propagator. Also, at finite m*, the 
static Coulomb field can be screened by the electrons at large distances. 

As we will discuss below, the situation in QCD is very similar - as /x — > 0, the periodic structure of the 6 vacuum 
becomes irrelevant. However, in the physical world /x ~ Aqcd ~ 200 MeV, so the ghost (describing the tunneling 
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FIG. 1: Fig. [T}a : Eq. © in momentum representation; helix lines represent gluons and the dashed line depicts the ghost. 
Fig. |T]-b : the gluon dressed by the interactions with the ghost: a “glost”, see Eq. (JTTJ for the corresponding self-energy 
expression. 


in the periodic 0-vacuum) strongly affects propagation of gluons with frequencies w ~ /z. At large distances, the 
ghost also gives rise to the screening of color charge, leading to the freezing of the effective coupling in the infrared 
IR limit for pure gauge theory, or to the vanishing of the coupling in the IR in QCD with light quarks. Because 
p ~ Aqcd ~ 1 / Rconf is on the order of inverse confinement radius R CO nf , these phenomena describe confinement of 
gluons at distances R CO nf — 1 fm. 

The key observation of our paper is that Eq. © and Eq. © define an effective ghost-gluon-gluon vertex (q,p). 
Using this vertex, we can re-write the correlator Eq. © at small q 2 as follows (see Fig. ©a): 

1 /* 1 ^ 

(<?) = J d A P (9,p) p 2( q _ p} 2 T " (®>p) = (9) 

From Eq. © we find that 

M 4 

(q,p) (q,p) oc , for q < p. (10) 

The vertices (q,p) describe the excitation of the ghost by gluons, and affect the gluon propagation at small 
virtualities. 

Indeed, the gluon propagator is now the solution of the equation shown in Fig. [2] with E (p) given by 


1 




(p) = (2tI)4 i j dAqT »i q 'P) -if S „\1 r -'i^p) 


q 2 (q-py 


( 11 ) 


where 1/q 2 is the propagator of the ghost. In evaluating the integral of Eq. (11) we assume that q -C p since T^’s 
describe the non-perturbative effects concentrated at small momenta (long distances). Therefore, 


(P) - [ dq 2 q 2 r M (g,p) -^T„(q,p) = - 4 f 

(27r) 4 J o q 2 p 2 p 4 J 


A 2 P 

dq = -guv-*, 

pZ 


( 12 ) 


where we used Eq. (10); note that in Eq. (12) we made rotation to the pseudo-Euclidean space. We can now write down 
the Schwinger-Dyson equation for the gluon propagator [25] G^ (p) = g^ G (p) in terms of (p) = g E (p), see 

Fig.@: 


g(p) = 4 + 4 s ip) °ip) 

pZ pZ 


(13) 


with the solution 


Gip) = 


1 


P 2 ~ E (p) 


P i- p 2 


(14) 
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FIG. 2: The graphic form of the equation for the gluon propagator, Eq. (131. The full propagator is denoted by the bold helix 
line; thin helix stands for the perturbative gluon propagator Go(p), and the blob for self-energy £(p) given by Eq. (Ill and 
Eq. p|. 


The propagator (141 has remarkable properties. First, G (p) has no infrared singularities and no gluon pole in 
the physical region. Indeed, this propagator has only complex poles at p 2 = ±ip 2 . As a result, gluons cannot be 


observed as particles in detectors - in other words, they are confined. Second, the propagator of the type of Eq. (14) 


was proposed by Gribov (2j as a solution to the problem of gauge copies - multiple solutions to the gauge fixing 
condition, see ©HD! for reviews. Hence we can state that introducing the coupling to the ghost (and thus taking 
account of the periodicity of the 9- vacuum) solves the problem of Gribov copies and leads to the confinement of gluons. 
The dimensionful Gribov parameter acquires a well-defined meaning of topological susceptibility Xtop = ft 4 related to 
the rf mass by the Witten-Veneziano relation 1113; since p ~ Aqcd, confinement emerges at distances of about 1 
frn. Note that close to the deconfinement transition, the topological susceptibility vanishes reflecting the restoration 
of symmetry [TTJQ2], see [13] for a review. Since at p —»• 0 the gluon propagator becomes perturbative, the 

restoration of Ua(X) symmetry and deconfinement should occur at the same temperature as suggested by the lattice 
data HD; however close to T c the non-perturbative interactions induced by p yf 0 are important. 


In our approach, the propagator Eq. (14) results from the admixture of the ghost to the perturbative gluon (see 
Fig- 0b), with an amplitude defined by the topological susceptibility /z 4 . We thus propose the following name for the 
particle with propagator given by Eq. (14) that represents a coherent mixture of a gluon and a ghost - a “glost” [2jJ. 


Unlike the ghost, the “glost” can be produced in a physical process, but unlike the perturbative gluon, it is confined 
and can propagate only at short distances ~ p~ 


Let us now re-consider the asymptotic freedom of QCD DUCES] using the “glost” propagator Eq. ([14|. In our 


derivation we will compute the interaction energy of two heavy quarks in the Coulomb gauge maim that is free from 
the Faddeev-Popov ghosts, so we can avoid dealing with two different types of ghosts. The dominant contribution 
responsible for the asymptotic freedom stems from the diagram of Fig. [3]-a. The contribution of this digram in 
perturbative QCD takes the form 


n (a) = ig 2 C\ 


d 4 k' 


(2tt) 4 z (fc - fc') 2 ( k ' 0 2 - fc 


/ _ (fc-fc') 2 y 

V k 2 k ' 2 r 


(15) 


n is related to E by E = k 2 n. 


For the “glost” propagator (141, Eq. (15) changes and takes a different form: 

lj(a) f _ (fc — k ) 2 _ (fcp 2 - fc 2 ) / (fc • k ) \ 

9 2 J W* ((fc - fc ') 2 ) 2 + /Z 4 ((fc ' 2 - fc ' 2 )) 2 + p^ fc 2 fc ' 2 ) 

1 1_f 1 }ra*(' 


= 3 g 2 C 2 [ S^Re 


Re^ 


(27t) 4 z 1 (fc - fc') 2 -t-z/z 2 J 1 k'£ - fc' 2 + ip 2 


fc 2 fc ' 2 - (fc • fc') 


(16) 


In Eq. (16) we have four terms obtained by choosing different signs of ip 2 . To illustrate the procedure of calculation, 
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let us evaluate one of them: 


U[ a) =3 g 2 C 2 


d 4 k' 


k 2 k' 2 - (k ■ k') 


(2ir) 4 i ((fc - k') 2 + ig?) {k' 2 - k 12 + i /x 2 ) fc 2 fc' 2 ' 
Introducing Feynman parameters a\ + a 2 + c *3 = 1 we obtain (P = k' — a\k and as = g 2 /An): 


(17) 


r* 1 —ai 


3 

y 2 fe 2 


dai 


da 2 


d 4 P 


k 2 P 2 - (fc ■ P)' 


integrating over fcp 




2^2 


integrating over P 3as „ 2 

' ^T° 2 


da.2 i 

V&2 J 

1-1 da2_ 

3 \/“2 , 


f*l — cxi 


dai 


J (2n) 4 i (a 2 kQ — P 2 — k 2 a i(l — ai) + i\x 2 [a\ + a 2 ))" 
d 3 P 3 P 2 


(2n) 3 16 ^ p 2 _|_ p * ai (]_ _ Q , 1 ) _|_ ^ 2 ( a i _(_ a2 ^ 


5/2 


r*l — OL2 


da i < — - + In 


2L 


i/fc 2 ai(l — ai) + ifj, 2 (ai + a 2 ) 


(18) 


where L is an ultraviolet cutoff in the integration over momentum. The integrals over a\ and a 2 can be taken 
analytically. The main features of Eq. (18) are the following: it has a logarithmic divergence at large k and is finite 
at k = 0. 

Summing all terms we find that we need to replace In {L 2 /k 2 ) of perturbative QCD in the diagram of Fig. §}a by 
the following function: 


In [L 2 /k 2 ) —» In (L,k,n) = 

f 16 


In 


8 J 

4 L 2 


da i 

\fa\ . 


pi — ai 


da 9 


+ In 


4L 2 


\jk 4 a 2 ( 1 - a 2 ) 2 + /x 4 (ai + a 2 ) 2 / \ \Jk 4 a 2 { 1 - a 2 ) 2 + /r 4 (ai - a 2 ) 


(19) 


It turns out that the same substitution has to be done in the expression for Fig. §b that in perturbative approach 
gives a positive contribution to the /3-function. However the contribution of the quark loop (Fig. |3|c) remains the 
same as in perturbative QCD. The sign of this contribution to the /3-function is also positive, and it leads to Landau 
pole and a “Moscow zero” [T8] . As a result the QCD coupling in our approach tends to zero in the infrared region 
of k —3- 0. 


If we choose the renormalization point k = fi, the running coupling takes the form 

as (m) 


a S (k 2 ) = 


1 + a s (/x) 


IliVc 


(In (L, k, g) — In (L, /z, /i)) - ^ In (k 2 /jj 2 ) 


( 20 ) 


In Fig. [4] we plot the coupling as as a function of k for two cases: QCD and pure gluodynamics in our approach 
(Fig. [4^a), and the comparison of as in gluodynamics with perturbative QCD calculations in the leading order (Fig.jdj- 
b). In both cases we choose the renormalization mass to be equal to the mass of the Z-boson; we use the value of 
fj, = 0.18 GeV from the original paper [4] where it has been determined from the mass of g'-meson, and N c = Nf = 3. 


One can see that replacing the gluon propagator by the propagator of the glost in gluodynamics removes the Landau 
pole and leads to the finite value of as at k = 0. On the other hand, with the inclusion of quarks, the QCD coupling 
vanishes at k = 0. At short distances, the running coupling is dominated by the perturbative contribution and so 
is not modified. The ghost affects the running coupling in a way that is quite different from the effect of a single 
instanton, which has been shown to increase the effective coupling at distances on the order of the instanton size 
113 EH- This may not be surprising as the ghost describes the effect of many instanton transitions throughout the 
0-vacuum. The screening effect of the ghost admixture is clearly a consequence of the fact that it is a spin-zero 
pseudoscalar “particle”. 


It has been argued by Dokshitzer J5T] that the experimental data indicate that in the IR region the QCD coupling 
remains effectively small: 

1 f- 1 ’ 

ao = — / dkas{k) ~ 0.5 for fij = 2 GeV. (21) 

AH J 
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In our approach we get ag 


0.59 for renormalization point k 


Mz, in reasonable agreement with Eq. (211. 


It is of fundamental interest to establish the microscopic dynamics responsible for the long-range correlations of 
topological current captured by the ghost. A recent study within the “deformed QCD” attributes these correlations 
to the topological order in the vacuum I22| . Our result suggests a link between the confinement and the long-range 
topological correlations in the QCD vacuum, and provides a practical way of computing power-suppressed corrections 
to QCD amplitudes, in particular the ones that are forbidden in the perturbative approach. 


To summarize, we propose to modify the gluon propagator in perturbative QCD by taking account of the periodic 
structure of the QCD 9- vacuum. Our prescription for the gluon propagator leads to the coupling of the gluons to the 
ghost saturating the anomalous Ward identity for topological current. The resulting “glost” propagator appears to 
have the functional form originally proposed by Gribov, in which the role of dimensionful parameter is played by the 
topological susceptibility \t 0 p = /x 4 . Our approach thus removes the Gribov copies that usually plague perturbation 
theory, and describes confinement of gluons at distances ~ /x _1 ~ 1 fm. We also find that the running coupling in 
the IR freezes in pure gauge theory, or tends to zero in QCD with light quarks. Because the topological susceptibility 
vanishes above the deconfinement transition, the “glosts” become usual perturbative gluons in the deconfined phase 
at high temperatures. The glost propagator leads to the exponential fall-off of the high-energy hadron scattering 
amplitude at large impact parameters needed to satisfy the Froissart bound; this can solve the long-standing problem 
of the perturbative approach in describing high energy scattering [2.'ll . In QCD amplitudes the coupling to the ghost 
can give rise to spin asymmetries |24j that are different from the usual perturbative approach - it will be interesting 
to study the resulting implications for spin physics at colliders. 
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FIG. 3: The first order corrections to the Coulomb energy. Helix lines denote the transversely polarized gluon, the double helix 
lines show the longitudinally polarized gluon, and the solid arrow lines depict the quarks. 
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a S ( k ) 




FIG. 4: The running constant as as a function of momentum k. Fig. [4 a shows our result for as in QCD with light quarks 
(blue line that goes to zero at k —» 0) and gluodynamics (red line). Fig. 4jb shows the comparison of momentum dependence 
of as in gluodynamics in our approach (blue line) and perturbative QCD (red line that is above the blue one at small k). The 
renormalization point is chosen at the mass of Z-boson, k = Mz■ 


















